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ABSTRACT 

The recent discoveries of X-ray and optical counterparts for GRBs, and a 
possible discovery of a host galaxy, implies that a direct measurement of the 
red-shift of some GRBs host galaxies is eminent. We discuss the implications of 
such measurements. They could enable us to determine the GRBs luminosity 
distribution, the variation of the rate of GRBs with cosmic time, and even, 
under favorable circumstances, to estimate G. Using GRB970508 alone, 
and assuming standard candles, we constrain the intrinsic GRB evolution to 
p{z) = (1 + 2;)-0-5±0.7_ 
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1. Introduction 


The recent observations of the Italian-Dutch Beppo/SAX satellite of 7 -ray 
bursts, (ICosta et al. 199'^) , with error boxes of a few arc-minutes across, enabled 
a follow-up by optical and radio observations and the discovery of x-ray, optical 
and radio counterparts to GRBs. The optical observations provided, for the hrst 
time, independent estimates of the distances to GRB sources, using absorption 
lines or association with host galaxies, ([Metzger, et al. 1997| ), and demonstrated 
beyond doubt the cosmological origin of GRBs. It is highly possible that in the 
months to come several GRBs would have independent red-shifts estimates. We 
may use these to obtain estimates of the luminosity diversity and the intrinsic 
evolution of the GRB rate with cosmic time. Under favorable conditions we 
might even be able to use GRBs as cosmic probes for estimating the cosmological 
parameters. 

Given a group of GRBs with measured fluxes and red-shifts we should 
first calculate the luminosity distribution of the sources . Then using this 
luminosity function we should estimate the theoretical peak-flux distribution 
in cosmological models. Gomparison of this distribution with the observed 
one would yield a direct estimate of the cosmic evolution of GRBs - that is 
the variability of the GRB rate per unit comoving volume per unit comoving 
time. Depending on the width of the luminosity function, the null assumption 
of no cosmic evolution could be proven or ruled out with a few dozen bursts. 
Gosmological parameters like the closure parameter, D, and the cosmological 
constant. A, influence only weakly the peak-flux distribution (jGohen fc Piranj 


1995|) hence this analysis can be done safely assuming D = 1 and A = 0 . 
However, if the luminosity function is narrow enough we might use the peak-flux 
red-shift relation to obtain a direct measure of D. 
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2. The Luminosity Function. 

The association of a host galaxy with a GRB provides us with the red-shift of 
the host, which is the bursts’ red-shift assuming that the burst is close to the host. 
Additionally we have the usual peak photon flux parameter, p{ph/cm?/sec), 
which is transferred later to apparent luminosity, l{erg/cm?/sec), using bursts’ 
spectra. The peak photon flux vs. red-shift and source luminosity is: 

^ N[ui{l + z)..M2{1 + z)] 

where vi = SOkeV, 1/2 = 300keV are the detector boundaries, 2 ; is the burst 
red-shift, c/H q is Hubble distance, and N[i'i,P 2 ] is the number of photons 
emitted in the range [i^i,r' 2 ]- We have used H = 1 and A = 0 in eq. |I|. 

The effect of H and A on the luminosity is not large, and we will discuss 
it later. The luminosity dependence on Hubble constant is a scale factor of 

When comparing bursts from different red-shifts one must recall that the 
observed peak-flux is in the range 50 — 300keV. This corresponds to different 
energy ranges at the sources. In order to discuss a single luminosity that 
classihes the burst, we will consider L = /so^v^ Lyi>di> at the source. To convert 
from the observed peak flux to the intrinsic luminosity we assume that the source 
spectral form is a power-law Ly = L(z//50keV)“"(2 — a)/(6 ^““ — l)/(50keV)^ 
in the energy range 50keV < u < 300keV(l -|- Zmax), so that wherever the 
source is, the detector sees a power-law spectra. We will use a = 1.5 for all 
bursts. This value is probably a good typical estimate ( [Band, et al. 199^ ), 
even though the spectra is not the same for all bursts . If necessary, one can 
use the measured spectrum of each burst to estimate its intrinsic luminosity at 
the 50 — 300keVband, however at this stage this simple estimate is sufficient. 
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Alternatively one can view the variability in the spectral index as an additional 
random variable that simply widens the luminosity function. Using this spectral 
shape we obtain: 


I = 


L{1 + z) 


L{l + z)-^hl 


75 


( 2 ) 


47r(2c/i/o)2(l - 7.7 * 105^(1 - 

We would be able to obtain a direct estimate of the GRB luminosity function 
in the observed band when as few as a dozen bursts will be observed. Using eq. 
^ we hnd the luminosity for each burst. Using those luminosities we can estimate 
the luminosity distribution function using maximum-likelihood or any other 
statistical method. Recall that current data, and in particular the peak-flux 
statistics of GRBs does not constrain the luminosity distribution of GRBs ( 
Gohen &: Piran 1995| , [Loredo fc Wasserman 19^ ). Glearly, if the sub-group 


of GRBs with optical counterparts is biased, this distribution estimate will be 
biased. 


3. Cosmological GRB Evolutions 


One of the interesting features that might distinguish between different 
cosmological GRB models is the rate that GRBs occur per unit time per 
unit comoving volume: p{z). These new measurements could yield a direct 
estimate of this distribution. Once the GRB luminosity distribution is known 
we can proceed and compare the theoretical peak-flux statistics (using the 
observed luminosity distribution) with the observed one. This distribution 
depends strongly on the intrinsic evolution of GRBs, that is on variation of 
p{z). Following ( Gohen &: Piran 1995|) we characterize this dependence as 
p{z) = (1 + z)~^. Gomparison of the theoretical and observed distribution would 
limit (3. The cosmological parameters G and A influence this distribution rather 
















weakly ( Cohen fc Piran 1995|) and have no substantial effect in estimating (3. 
Recall that prior to the independent knowledge of GRBs luminosity, one could 
not distinguish between cosmological effects and intrinsic evolution using count 
statistics. 

In fact this comparison can be done even with the current data and assuming 
a narrow luminosity distribution (standard candles). We can use GRB970508 
to constrain the evolution. Using the peak flux = 1.6 * ergjsecjcw?^ 

( IKouveliotou, et al. 1997|) , the red-shift of the absorption lines 2 ; = 0.835 
( pVIetzger, et al. 1997|) which sets a lower limit 2 ; > 0.835 for the burst, and 
the absence of prominent Lyman-alpha forest in the spectrum which compose 
an upper limit 2 ; < 2.1, we obtain (3 = —0.1 ± 1.3 in 99% confidence level. 
Assuming that the absorption line of GRB970508 correspond to its own red-shift 
we estimate /3 = 0.5 ± 0.7 with this confidence level , see fig 1. The simplest 
hypothesis of no evolution /? = 0 is consistent with the observations. A milder 
assumption of Gaussian luminosity distribution with ai = Lobs/‘^, instead of 
standard candles, gives a lower limit (3 > —0.7 and no upper limit. 


4. Estimates of Cosmological Parameters. 

Despite numerous attempts to estimate the cosmological parameters, there 
are still large uncertainties. One may wonder whether GRBs would provide 
a meaningful independent estimate of these parameters. Using GRBs count 
statistics alone, D could not be estimated from the current data (|Cohen fc Piran 


1995). However, given a cosmological distribution of sources with measured 


red-shifts, we can try to estimate the cosmological closure parameter, D, in 
a similar manner to the attempts to estimate D from type I supernovae by 
Perlmutter, et al. 1996| . 
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The observed peak-flux depends on hi as: 


/ = /(L, (]) 


L 


__ 

UT^{zn/2 + {yt/2 - i)(v/nJTi -1))2 


( 3 ) 


Using the known parameters of each burst (peak-flux and red-shift) we obtain 
for each burst a function Lj = Li{VL). In the previous sections we have assumed 
n = 1 and obtained Lj(r2 = 1) ). For standard candles all Lj must be equal. 
Given two sources we have two functions L = Li^ 2 (G), with two variable and we 
should be able to determine G. 

A luminosity distribution will induce an uncertainly in this estimate that 
can be approximated by: 


dVt 


^ ^_ (Vl + z — l)Vl -|- 2 : _^ 

2 ^/2 — 3/2-\/l -|- 2 ' -7 3/2 -|- -|- z) L 


( 4 ) 


For 2 : = 1.5 we obtain ctq ~ 2(Tl/L. Thus assuming ul/L = 1, we need 100 
bursts with a measured 2 : to estimate G with an accuracy of = 0.2. Such a 
goal could be achieved within several years. At present it is not known whether 
the GRB luminosity distribution is narrow enough and satisfies this condition. 
However, as we have shown at section 2, the width of the luminosity distribution 
will be known soon. 


5. Future Detectors 

In view of the promising avenues that these observations have opened it is 
worthwhile to examine what will be the effect of future more sensitive detectors 
on these estimates. Gohen and Piran (1995) have estimates that BATSE is 
sensitive to bursts up to Zmax = 2. A more sensitive detector ( by a factor of 10 
) will detect bursts up to 2 : = 6.9. Assuming that p{z) is constant up to high 
value of 2 we And that the number of observed bursts will increase only by a 



factor of 2.1. The results are slightly better if the current z^ax is smaller. For 
example a ten fold more sensitive detector will measure 2.6 times more bursts 
than BATSE if Zmax=^-'^ and 3.5 times more bursts if Zmax = 1- 

At hrst sight these results might look discouraging. However, the rate of 
GRBs at high red-shift is unknown and could be critical in determining the 
nature of GRBs. At present it is not known if there are bursts which originate 
from high 2 ;. Most models that are based on compact objects cannot produce 
sources at very early time. On the other hand these model, and in particular 
the neutron star merger model predict a high rate of GRBs that will follow with 
a time lag of up to 10® years any extended star formation activity. This is 
approximately the time it will take the stars in a binary system to hnish their 
life cycle, become NS, loose angular momentum through gravitational radiation, 
and merge. Nuclear abundances measurement indicates that heavy elements 
that are produced in supernovae began to be produced earlier than heavy 
elements produced in neutron star mergers. It will be intriguing to see whether 
GRB rates would follow the trend of supernovae or neutron star mergers. 

Furthermore, , one has to recall that the relevant question for our purposes 
is not how many bursts are observed but how many bursts are observed with 
measured red-shifts. Gurrently the rate of detection of bursts with counter-parts 
is about one per month and from those detected until now only one has a 
measured red-shift. Here there is an enormous potential for improvements. 

For example systematic measurements of the red-shift of all bursts observed 
by BATSE 300 per year) would yield an independent estimate of G, with 
(Tq = 0.1, even if the luminosity function is wide, {ul/L = 0.9), within one year. 
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6. Discussion and Conclusions 


As expected, the direct red-shift measure of GRB970508 agrees well with 
estimates made previously using peak-flux count statistics (|Fenimore, et al 


19931, [Loredo fc Wasserman 199^ , [Cohen fc Piran 1995|) . It is remarkable what 


could be done with even several additional red-shifts. An estimate of the bursts’ 
luminosity distribution can be obtained even with a few bursts. This luminosity 
function combined with the observed peak-flux distribution would provide us 
immediately with an estimate of the cosmological evolution of the rate of GRBs. 

It is generally accepted that Fireball is inevitable in any cosmological model. 
Within this model the observed y-rays are produced during the conversion of a 
relativistic energy flow to radiation. However, the source itself that produces the 
flow remains unseen ([Piran 1996|) . The limits on cosmological evolution could 
shed light on the GRB mystery, by distinguishing between different cosmological 
models. For example the expected rate of merging neutron stars depends on the 
red-shift in a drastically different way than the expected rate of events related 
to AGNs. 

The implications of these red-shift measurements could reach even further. 
Detection of a signihcant group of GRBs with red-shifts could enable us to 
utilize GRBs to study cosmology. Hundred bursts with associated red-shifts will 
enable us to estimate the cosmological parameter G even if the GRB luminosity 
distribution is relatively wide. This goal is not practical with the current 
detection rate. However, if the luminosity distribution is narrow, or if a novel 
detection technique could be found which will yield a signihcantly higher rate 
of counter-part detection, we would be able to measure G using this method 
within several years. 
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Fig. 1.— The likelihood function (levels 31.6%, 10%, 3.16%, 1%, etc. of the maximum) in 
(/3, L) plane for standard candles., a = 1.5, = 1, and evolution given by p{z) = (1 + z)~^. 

Superimposed on this map is the luminosity of GRB970508, with thick line where the 
likelihood function >1%. We have used = 1. 




